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A novel approach to identifying 
antiarrhythmic drug targets 



Sudden cardiac death, secondary to ventricular fibrillation (VF), remains 
the leading cause of death in the USA. Recent experimental and theoretical 
studies suggest that VF could be caused by spiral wave re-entry. The 
initiation and subsequent break-up of spiral waves has been linked to 
electrical altemans, a phenomenon produced in cardiac tissue that has a 
steeply sloped restitution relation. Agents that reduce the slope of the 
restitution relation have been shown to suppress altemans and, presumably 
by that mechanism, terminate VF. These results suggest that electrical 
restitution could be a promising new target for antiarrhythmic therapies. 

▼ Catastrophic rhythm disturbances of the 
heart, such as ventricular fibrillation (VF), are 
the major cause of death In the United States 
[1,2]. Treatment of these rhythm disorders re- 
mains largely empirical, in part because of an 
incomplete understanding of underlying cel- 
lular electrophysiological mechanisms. Recent 
studies, building on earlier theoretical work by 
Krinsky, Wlnfree and colleagues [3,4] and exper- 
iments by AUessie [5], have suggested that spiral 
wave re-entry could be the 'engine' that drives 
VF [6-11]. Although there Is substantial evidence 
that spiral wave re-entry contributes signifi- 
cantly to the induction and maintenance of 
VF, the exact mechanism by which spiral waves 
sustain VF is currently being debated. Correctly 
identifying the mechanism (or, more likely, 
mechanisms) by which spiral waves cause VF is 
likely to be a crucial step in the development of 
pharmacological approaches to VF prevention. 

Restitution hypothesis for VF 

One of the hypotheses to account for the ap- 
parent link between spiral waves and VF is the 
restitution hypothesis, which proposes that 
VF is caused by the break up of a single spiral 
wave into multiple self-perpetuating wavelets 
[8,10,12-17]. This mechanism is similar to 
that proposed decades ago by Moe [18], but 
with one important difference: it can occur in 



intrinsically homogeneous cardiac tissue (for 
reasons to be described below). The transition 
from normal planar wave excitation to a single 
spiral wave and, ultimately, to multiple wavelets 
is thought to underlie the transition from 
normal sinus rhythm to ventricular tachycar- 
dia and VF characteristic of patients who suc- 
cumb to sudden death (Fig. 1). 

The exact mechanism for the break up of 
spiral waves is unknown. However, there is 
considerable evidence that break up is closely 
related to action potential duration (APD) resti- 
tution, which is the relationship between APD 
and diastolic interval (Dl, the time interval be- 
tween action potentials), in which APD is de- 
termined by the preceding DI (Fig. 2). During 
pacing of cardiac tissue at progressively shorter 
cycle lengths, APD decreases, until at suffi- 
ciently short cycle lengths a period-doubling 
bifurcation occurs and APD begins to alternate 
between a long duration and a short duration, 
a phenomenon known as APD altemans (Fig. 2) 
[19-24]. The development of APD altemans 
requires that the slope of the restitution relation 
exceeds 1, whereas if the slope of the restitution . 
is less than 1, altemans will not occur. 

The slope of the restitution relation, and the 
corresponding presence or absence of APD alter- 
nans, has been linked to spiral wave stability. If 
the slope of the APD restitution relation is <1, a 
spiral wave tends to stabilize and produce a peri- 
odic rhythm, the manifestation of which might 
be monomorphic ventricular tachycardia. If, 
however, the slope of the APD restitution rela- 
tion is £1, a single spiral wave might disintegrate 
into many spiral waves, manifested as VF. These 
transitions are illustrated in Fig. 3, using results 
generated by a computer model Once initiated, 
a single stable spiral wave can be destabilized by 
increasing the slope of the restitution relation 
from <1 to H. Similarly, multiple wavelets can 
be induced to coalesce into a single spiral wave 
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by reducing the slope of the restitution 
relation from Si to <1. Thus, a reduction 
f the slope of the restitution relation 
from Si to <1 is expected to prevent 
the induction of VF and to convert 
existing VF into a periodic rhythm. 

Several recent studies have provided 
experimental and theoretical support 
for a causal relationship between APD 
restitution and VF [12,13,25-29]. If the 
APD restitution relation contains a 
region of slope Si, APD altemans and 
VF are induced by pacing at short cycle 
lengths. If the slope of the APD restitu- 
tion relation is reduced to less than 
one, APD alternans is suppressed and 
VF is not induced. Furthermore, if a 
fibrillating ventricle is exposed to a 
drug that reduces the slope of the resti- 
tution to less than 1, VF Is converted 
to a periodic rhythm sustained by a 
single stable spiral wave [12,13,26]. 

Unfortunately, interventions that sup- 
press VF experimentally (verapamil [13], 
bretyllum [12] and hyperkalemia [26]) 
have unwanted actions that severely 
limit their clinical use. Nevertheless, 
the effects of these interventions on 
restitution and VF have provided valu- 
able insights regarding new, potentially ' 
more clinically relevant, drug targets. 
In particular, attempts to understand 
the mechanism by which calcium 
channel blockers alter APD restitution 
have led to the realization that increas- 
ing selected outward repolarizing cur- 
rents might also flatten restitution, as 
discussed in the following sections. 

Caldum channel antagonists and VF 

The suppression of VF by verapamil, 
in association with a reduction in the 
slope of the APD restitution relation 
(Fig. 4), suggests that the i-type caldum 
current (Iq) has a key role in restitution 
and in the development of VF [13,30]. 
To better define these roles, we recently 
developed an ionic model for the cardiac 
ventricular myocyte (CVM) based on the 
Window [31] and Luo-Rudy models [32] 
that generates physiologically realistic 
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Figure 1 . Proposed evolution of cardiac wave propagation patterns underlying the 
transition from a normal cardiac rhythm to ventricular fibrillation. Transition from normal 
sinus rhythm (NSR) to ventricular tachycardia (VT) and ventricular fibrillation (VF), as 
recorded on the surface electrocardiogram. Corresponding transitions from a planar 
wave to a single spiral wave to multiple spiral waves. 
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i of action potential duration (APD). (a) Definition of APD, 
diastolic interval (Dl) and cycle length (CI). The duration of the n+1th action potential 
(APD w1 ) is a function of the preceding Dl (Dl„). (b) Plot of APD versus preceding Dl. 
(c) Plot of APD versus CL Crosses are APDs during the CLs in which no altemans 
occurs. During APD alternans, unfilled circles are APDs during the long action 
potential, and filled circles are APDs during the short action potential. Modified, with 
permission, from Roller et al. [25]. 





Figure 3. Wave patterns in a computer model of cardiac activation. Sequential snapshots 
(A-)) of activation during the development of a stable single spiral wave under conditions 
in which the slope of the restitution relation for action potential duration is <1 (left panel). 
The transition from a single spiral wave into many spiral waves for a restitution slope 
>1 (A-O, with some intervening steps deleted: right panel). Results courtesy of 
Dante R. Chialvo (Department of Physiology, University of California, Los Angeles: 
http://www.uda.edu). 




Figure 4. Effects of verapamil on action potential duration (APD) restitution and 
ventricular fibrillation (VT). (a) Example of the reduction in the slope of the APD 
restitution relation after exposure to verapamil, (b) Corresponding reduction in the 
magnitude of APD altemans (AM) during pacing at different basic cycle lengths (BCL) 
after exposure to verapamil, (e) Effects of verapamil on the composite Fast Fourier 
Transform (FFT) of monophasic action potential (MAP) recordings during VF in arterialry 
perfused canine ventricle. Verapamil, added after 35 min of VF had elapsed, converted 
the broad FFT spectra during VF to a single frequency. Modified, with permission, from 
Riccioeto).[13]. 



APD altemans over a wide range of 
pacing basic cycle lengths (BCL) [33]. 
Using this model the L-type Ca 2 * current 
was implicated as an important determi- 
nant of APD altemans, according to the 
mechanism shown in Fig. S. On initi- 
ation of pacing at a short BCL, Iq, is fully 
recovered before the first action poten- 
tial and activates fully during the 
action potential, which results in a long 
APD. The long action potential is subse- 
quently followed by a short diastolic 
interval, during which Iq, fails to recover 
completely from Ca^-induced inacti- 
vation. Because of decreased availability 
of Ic the duration of the subsequent 
action potential is shorter, which results 
in a longer succeeding diastolic interval, 
more complete recovery of 1^ and a long 
action potential duration. This cycle 
then repeats, eventually establishing a 
steady-state altemans of Ic, and APD. 

Potassium channel agonists and VF 

From the mechanism for APD altemans 
outlined in Fig. S, we anticipate that 
shortening of APD during rapid pacing 
will prolong DI, which, in turn, might 
provide adequate time for complete 
recovery of Ic,. If so, reducing the 
magnitude of APD altemans could be 
accomplished by increasing outward 
repolarizing currents, rather than by 
decreasing Lo,. This idea has been sup- 
ported by additional computer model- 
ing studies in which APD altemans 
was suppressed by increasing any one 
of several outward repolarizing cur- 
rents [e.g. the rapid (I*,) and slow (Ik») 
components of the delayed rectifier 
and the inward rectifier (I K ,)] {33]. 
Although agonists for these currents 
are generally not available, it might be 
possible to modify currents such as Ik, 
and Ik, by increasing phosphatidyl in- 
ositol bisphosphate (PIP2) levels [34] 
or by altering the phosphorylation state 
of the channels [35-37], provided the 
channels can be phosphorylated with- 
out concomitant phosphorylation of 
calcium channels or can be upregulated 
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Figure 5. Putative mechanism for action potential duration (APD) 
altemans and the suppression of attemans by calcium channel 
blockade, (a) Under control conditions, alternation of the 
length of the diastolic intervals (Dl) during APD altemans causes 
alternation of the magnitude of Ic^ secondary to Incomplete 
recovery of l a from inactivation during the short Dl. (b) Reduction 
of Ic decreases APD, which prolongs Dl, so that Dl following 
each APD is sufficiently long to enable complete recovery of 1^ 
In the absence of alternations of 1^ APD altemans is suppressed. 



in the presence of a calcium channel antagonist, to offset 
increased 1^ secondary to phosphorylation [38]. 

Recently, we have focused on increasing I fc as an ap- 
proach to suppressing APD altemans. 1^ has an important 
role in cardiac repolarization, increasing to a maximum 
during phase three of the action potential, as the channel 
recovers from inactivation, and then decreasing as the 
electrical driving force decreases and as deactivation of the 
channel increases [39-43]. Because I fc contributes mini- 
mally to the action potential plateau, increasing I & can 
have little or no effect on the Ca 2 * transient. Consequently, 
by increasing Ir, it can be possible to suppress APD alter- 
nans without adversely affecting contractility. 

At present, however, there are no I fo agonists available 
to test this hypothesis. To circumvent this problem, we in- 
fected isolated ventricular myocytes with an adenovirus 
expressing HERG, the gene that encodes the pore-forming 
domain of I*,, to increase HERG protein expression level 
and the corresponding l to current [44] [Hua, F. et al., 
unpublished data]. After verifying that I*, had increased 
post-infection, the myocytes were paced at rapid rates to 
determine whether increasing I*, by overexpressing HERG 
suppressed APD altemans. As shown in Fig. 6, overexpres- 
sion of HERG markedly increased I*,, as recorded during 
action potential clamp, and suppressed APD altemans 
during rapid pacing. Overexpression of HERG did not, 
however, significantly reduce I 0 , suggesting that this ap- 
proach to suppression of altemans need not be accompa- 
nied by a reduction in contractility. 

The current paradigm regarding I*, arid VF is that block- 
ing Ik, is expected to be anti-arrhythmic, secondary to 
prolongation of APD and refractoriness [45,46]. However, 
this approach has not been successful in preventing VF 
and, moreover, is associated with pro-arrhythmia [45]. 
Prolongation of APD, as reflected by an increased duration 
of the QT interval on the EGG, and an increase in the inci- 
dence of ventricular arrhythmias are also associated with 
the administration of numerous non-cardiac drugs that 
block Ir, [47-50]. Similarly, inherited loss-of-function 
mutations in l b are accompanied by a prolongation of the 
QT interval and by an increased risk of lethal ventricular 
tachyarrhythmias, such as torsade de pointes [51-54]. 

The cardiac arrhythmias associated with inherited or 
drug-induced abnormalities of 1^ are thought to be precip- 
itated primarily by bradycardia-induced prolongation 
of repolarization [51,52]. The potential mechanisms by 
which Ik, might facilitate the induction of arrhythmias at 
slow heart rates have been studied extensively [46,50,53]. 
However, the contribution of l & to repolarization during 
tachycardia, which might also be important for the 
development of cardiac arrhythmias, has not been well 



characterized. The observation that reducing I*, increases 
the magnitude of APD altemans could provide an addi- 
tional mechanism to account for the proarrhythmic effects 
of Ik, blockers, in that increased altemans magnitude 
would be expected to destabilize ventricular tachyarrhyth- 
mias, leading to the development of VF. Conversely, the 
observation that increasing I b reduces the magnitude of 
APD altemans and the slope of the APD restitution rela- 
tion provides a rationale for the development of a new 
class of compounds, 1^ agonists, with the expectation that 
such compounds might have anti-fibrillatory effects. 

The latter hypothesis relies, however, on the expectation 
that Ik, can be increased sufficiently to reduce the slope of 
the restitution relation without shortening action po- 
tential duration to such an extent that contractility is 
impaired, secondary to a shortening of the plateau dura- 
tion and attenuation of Ic,. In addition, shortening of APD 
could lead to a reduction in the wavelength of re-entry 
circuits (in which wavelength = refractory period x conduc- 
tion velocity). Reduction of the wavelength could, in turn, 
precipitate wavebreak by other mechanisms (e.g. so-called 
'head-tail' interactions, in which a wavefront encounters a 
waveback and fragments) [55]. These and other potential 
drawbacks to the use of K -channel agonists to suppress VF 
remain to be evaluated critically. 

Potential implications for drug development and 
evaluation 

Until recently, therapy for the prevention of sudden cardiac 
death had been based on the presumption that frequent 
ventricular ectopy, in particular ventricular tachycardia, is 
a prelude to ventricular fibrillation [2]. Accordingly, drugs 




Figure 6. Effects of HERC overexpression on I* and action potential duration (APD) 

i isolated canine cardiac myocytes, (a) Current traces of l„ (upper panel) and 



action potentials during pacing at CL = 200 msec (lower panel) obtained from myocytes 
in the day of isolation (Day 0). (b) Current traces of I* (upper panel) and action 
potentials during pacing at cyde length (CL ) = 200 msec (lower panel) obtained from 
myocytes after three days in cell culture and infection with adenovirus (Day three + HERC). 
HERC overexpression increased Ik, and suppressed APD altemans. 



that does not involve alteration of resti- 
tution kinetics (e.g. slowing of conduc- 
tion oi prolongation of refractoriness). 

Given these expectations, the value 
of existing agents for the prevention of 
VF could be re-evaluated in the light of 
their effects on the restitution relation 
and new drugs targeted against resti- 
tution could be developed. Naturally, 
such an effort would need to recognize 
that factors other than APD restitution 
significantly contribute to the induction 
and maintenance of VF. Furthermore, 
drugs intended to alter restitution se- 
lectively might have additional effects 
that could offset their intended effect. 
Nevertheless, judicious alteration of the 
APD restitution slope appears to be a. 
promising approach to the treatment 
of VF and, as such, represents a poten- 
tially fruitful opportunity for drug 
development. 



that suppress inducible or spontaneously cxxurring ventricu- 
lar tachycardia are expected to prevent sudden death. 
However, recent large-scale clinical trials have indicated 
that classes of drugs that are effective for the suppression 
of ventricular tachycardia, for example, Class I and Class III 
antiarrhythmic drugs, do not prevent sudden death and, 
in fart, could be proarrhythmic [45,56]. In contrast, other 
classes of drugs that are not particularly effective for the 
suppression of most forms of ventricular tachycardia, such as 
fj-adrenergic receptor antagonists [57] and calcium channel 
antagonists [58], could reduce mortality from sudden death. 

If a causal relationship between the slope of the APD 
restitution relation and VF is confirmed, such a result 
could have significant implications for the pharmacologi- 
cal therapy of sudden cardiac death. Drugs that reduce the 
slope of the restitution relation would be expected to pre- 
vent the development of VF, but would not be expected to 
suppress ventricular tachycardia, if ventricular tachycardia 
is caused by some variant of spiral wave re-entry. Conversely, 
drugs that do not reduce the slope of the restitution relation 
would not be expected to prevent VF, although they might 
suppress ventricular tachycardia, perhaps via a mechanism 



Conclusions 

Electrical restitution is a recent addi- 
tion to factors that play a key role in 
the development of ventricular tach- 
yarrhythmias. Experimental interven- 
tions that reduce the slope of the resti- 
tution relation have been effective in suppressing VF. These 
results encourage further investigation of altering restitu- 
tion as a means of preventing sudden death, with the real- 
ization that clinically useful interventions that flatten 
restitution without having untoward effects have yet to be • 
developed. 
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Dynamic mechanism for conduction block in 
heart tissue 



Abstract Previous work has shown that dynamic heterogeneity and conduction 
block can occur in homogeneous heart fibres during prolonged pacing at rapid 
rates. Here we investigated the mechanism for conduction block following the 
delivery of one to four premature stimuli using a coupled maps computer model 
of a one-dimensional canine heart fibre. The coupled maps model allowed 
us to identify the roles that velocity (V) restitution, action potential duration 
(D) restitution and cardiac memory (M) played in the development of spatial 
heterogeneity and conduction block We found that the likelihood of conduction 
block could be reduced by three methods. ( 1) By altering the V restitution function 
so that conduction slowed at very short rest intervals (/). (2) By altering the D 
restitution function to reduce the sensitivity of D to changes in / . (3) By increasing 
the contribution of cardiac memory (Af)- Although the results of this study need 
to be confirmed experimentally, they suggest several potential interventions that 
may reduce the probability of arrhymmia induction. 
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1. Introduction 

Catastrophic heart rhythm disorders are among the leading causes of death in the United States. 
The most dangerous of these arrhythmias is ventricular fibrillation, a disturbance in which 
disordered wave propagation causes a fatal disruption of the synchronous contraction of the 
ventricle. Although the exact mechanism for fibrillation is still being debated, one theory 
proposes that fibrillation is a state of spatiotemporal chaos consisting of the perpetual nucleation 
and disintegration of spiral waves [4,26], in association with a period doubling bifurcation of 
local electrical properties [13, 14, 16, 18]. Nucleation of the initiating spiral wave pair is caused 
by local conduction block (wave break) secondary to spatial heterogeneity of refractoriness in 
the ventricle [2, 4, 22, 25, 26]. Until recently, spatial heterogeneity was thought to result solely 
from regional variations of intrinsic cellular electrical properties [22, 27] or from stimulation at 
more than one spatial location [20, 24, 25]. However, it is now appreciated that purely dynamical 
heterogeneity can be sufficient to cause conduction block during single-site stimulation in both 
homogeneous one-dimensional models of canine heart tissue and in rapidly paced canine Purkinje 
fibres [8, 12], A similar mechanism has been shown to precipitate conduction block and spiral 
break-up in models of homogeneous two-dimensional tissue [9]. 

The period doubling bifurcation implicated in the transition to conduction block is 
manifest as altemans, a beat-to-beat long-short alternation in the duration of the cardiac action 
potential [13, 14, 16, 18,24]. Previous investigators have hypothesized that altemans can be 
accounted for by a simple um-dimensional return map called the action potential duration 
restitution function [5, 6, 15, 17]. This hypothesis assumes the duration D of an action potential 
depends only oh its preceding rest interval / through some function /(/) that is measured 
experimentally. If the D restitution function has a slope > 1, then a period doubling bifurcation 
occurs for some value of the stimulus period T, where T = D + /. The velocity V at which an 
action potential propagates can also be described by a restitution function, where V = c(7). 

It has been shown previously that the combination of a steeply sloped D restitution 
function and a monotonically increasing V restitution function is sufficient to produce dynamical 
conduction block during sustained pacing at a short cycle length [1 1, 12]. This observation 
may provide a generic mechanism for wave break and the onset of ventricular tachycardia and 
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fibrillation. However, it is unlikely that the conditions used to demonstrate this phenomenon 
experimentally apply to the clinical situation, where the induction of ventricular tachyarrhythmias 
is typically associated with the interruption of normal cardiac rhythm by only a few premature 
beats. A single premature beat is sufficient to cause spatial heterogeneity in the form of discordant 
altemans [24], but the conditions required for the development of conduction block m this setting 
have not been investigated. 

To address this issue, in the present study we determined whether dynamic heterogeneity 
and conduction block occur in a computer model of canine heart fibres in which pacing at a slow 
rate is interrupted by one to four premature stimuli. This protocol simulates the interruption 
of sinus rhythm by one to four premature ventricular complexes (PVCs), a situation that can 
lead to the onset of ventricular fibrillation clinically. Furthermore, we determined whether such 
conduction block can be accounted for by the same dynamical mechanism that underlies the 
development of conduction block during rapid pacing. A coupled maps model of heart tissue 
was used for the study, which allowed us to identify the roles that V and D restitution play in 
the development of dynamic heterogeneity. We also considered the additional contribution of 
. cardiac memory (M) to the development of dynamical heterogeneity and conduction block. 



We studied the mechanism for conduction block following the delivery of multiple premature 
stimuli using a one-dimensional computer model of a canine heart fibre. The length of the fibre 
was 4 cm, which is similar to the length of fibres used for previous experimental studies [12]. The 
fibre was stimulated at one end at a cycle length near normal canine sinus rhythm (S 1 = 500 ms). 
After ten beats at SI, four premature stimuli were delivered (S2, S3, S4 and S5) at the same site. 
The coupling intervals between these premature intervals were varied and the conduction of the 
resultant action potentials was observed (figure 1). For each combination of premature stimuli, 
one of three possible outcomes occurred: 

(1) the stimulus elicited an action potential that propagated down the entire fibre; 

(2) the stimulus did not produce an action potential (type I block); 

(3) the stimulus elicited an action potential that blocked before reaching the end of the fibre 
(type H block). 

The first two cases would not be conducive to die development of wave break and initiation 
of re-entry, in that conduction of the premature response either does not occur at all or occurs 
equally well everywhere along the fibre. The third case, however, could lead to wave break 
and spiral wave initiation if it occurred in two- and three-dimensional tissue, provided the block 
was local. The development of local block would be facilitated in intact myocardium by twist 
anisotropy and intrinsic heterogeneity. 

To assess the vulnerability of the simulated tissue to type II block, the S2-S3-S4-S5 
combinations were applied in the following way: the S1-S2 interval was varied from the 
nunimum value that conducted (Sl-S2 min ) to Sl-S^ + 20 ms. For each S1-S2 interval, the 
S2-S3 interval and the S3-S4 interval were varied in combination from the minimum value that 
conducted up to a value of 250 ms for each interval. Finally, foreach S2-S3-S4 combination that 
conducted, the S4-S5 interval was varied from the minimum that generated an action potential 
at the site of stimulation to the minimum value that conducted down the entire fibre. If no 
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Figure 1. (a) Schematic representation of the relationships between stimulus 
period (T), actioa potential duration (D), diastolic interval (/) and memory (Af). 
(b) Example of the stimulation protocol. The model fibre was paced for ten beats 
at a constant Sl-Sl interval, after which a series of premature stimuli (S2-S5) 
was delivered. R indicates the refractory period of a cardiac cell, during which 
the delivery of a stimulus does not produce a propagated response. 

S2-S3-S4-S5 combination was found that produced conduction block for more than a 20 ms 
window in either the S2-S3 interval or the S3-S4 interval, the search was halted for that interval. 
All intervals were varied in steps of 1 ms. 

2.1. Coupled maps model 

The study was conducted using a coupled maps model of a one-dimensional cardiac fibre that 
has been described in detail elsewhere [10, 12]. Briefly, the model is based on the equation 

W*) = T^i(xi) - A. + i(x,). (1) 
r n+1 (*, ) was the time interval between activations of site x, . It was determined by including the 
time delays caused by the propagation from the pacing site to site x t , which yielded 

£^ Ax Q Ax 0 . 

r was the time interval between activations applied to the pacing site and Ax = 0.1 was the 
length of a single cell (time units in milliseconds and space units in millimetres). The conduction 
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velocity V B (x,) depended only on / through the velocity recovery function V„ = c(/„) given by 
c(J) = Wl - exp(-(/ + ?)/&)). = 0.72, 5 = 14, and was, varied to adjust the value 
of V at / = /nun- & was determined locally based on a memory model mapping [10] given by 

Af„ + , = g{M n , /„, A,) = e-^ T -[l + Wn ~ l)t~ D ' ,X °) 

/ B \ O) 

D n+1 = /(AU, /„) = (1 — ciM n +\)\A + — fcToT^ j- 

was the time constant of accumulation and dissipation of memory (both constants were chosen 
to be the same). . A = 88, B = 122, C = 40 and r m = 180 [10]. o and x D were varied in this 
study, a (which varied between zero and unity) determined the influence of memory on D, and 
x D was used to adjust the dependence of the D recovery function on /. We note that, for memory 
models, the dependence of D on / (as well as the occurrence of a period doubling bifurcation) 
is not related to the steady-state restitution slope in a simple manner [10]. 

Coupling between sites was included by using the diffusion terms from Echebarria and 
Karma [8]. These terms modelled the electrotonic current that flowed out of (into) a cell and 
into (out of) its neighbour if the action potential of the first cell was longer (shorter) than that of 
its neighbour. Including diffusion then yielded 

D n+l = /(Mn+i, /„) + ^D^x - wVD n +u (4) 

with £ = 1.0 and w = 0.35. Discretizing the derivatives produced a tri-diagonal linear system 
of equations thatcould then be solved easily. The defining equation for the model was therefore 

• M B+ , = g(M„, /„, D n ) 

Z) B+ i=/(M B+1 ,/ B ) + ? 2 V 2 Z>„ +1 -u;VZ) B+1 

Finally, conduction block was modelled by setting / = 0 for / < 1^ = 2. 
3. Results 

3.1. Characteristics of conduction block in the wild type model 

Figure 2 shows a typical example of type II conduction block in the baseline model (hereafter 
called the wild type model). Despite the homogeneous state that existed at the end of the SI 
stimuli, dynamical heterogeneity developed following delivery of S2, secondary to V restitution 
(see below). Heterogeneity was manifest as a short-to-long gradient in / , which was magnified by 
subsequent premature stimuli until the action potential generated by the S5 stimulus encountered 
a region where /</„„, and conduction failed. 

Figure 3 shows a histogram of type II conduction block at various S2 and S3 intervals 
generated from the wild type model. Note the presence of a large peak centred at an S2 very 
close to the minimum value for conduction and at an S3 approximately 50 ms longer than the 
minimum value for conduction. Surrounding the peak is a flat plateau region of five to ten blocks 
per bin that extends throughout the entire (S2, S3) region that was explored numerically. An 
ovemhelming majority of these instances of block occurred after the S5 stimulus. However, 
six examples of 'early block' were found after an S4 stimulus, and one was found after an S3 
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" Sl=500ms S2=207ms S3=146ms S4=148 ms S5=106ms 

Figure 2. Example of type II conduction block in the coupled maps model. Each 
panel shows / as a function of cell number for the last SI = 500 ms beat Geft 
panel) and for beats S2 = 207 ms, S3 = 146 ms, S4 = 148 ms and S5 = 106 ms. 
Note the difference in scale for each beat Conduction block, corresponding to 
/ < /„nn, occurs near cell 100 on the S5 beat and is denoted by two horizontal 
black lines. 

stimulus. The latter explain the presence of the plateau in the histogram; the plateau is made up of 
'degenerate' counts in which long S2 or S3 intervals that maintained the initial SI homogeneity 
were followed by early block intervals leading to block at S5. 

Since conduction block occurs when / < changes in 7^ will produce significant 
changes in the incidence of conduction block. Because is a component of the functions that 
control V restitution, D restitution and M, changes in U influence all three dynamic variables. 
In this study, we held U fixed and individually studied the roles that V restitution, D restitution 
and M played in detennining the risk of type H block, as assessed by the S2-S3-S4-S5 protocol. 

3.2. Contribution of V restitution to conduction b lock 

V restitution played two important roles in the development of conduction block First, as 
mentioned above, V restitution was the source of the initial heterogeneity after delivery of the 



101.7 




Figure 3. Histogram of the incidence of type II conduction block in the wild type 
model, a - 0.2, t D = 28 and p = 17.408. The height at each point in (S2, 
S3) space corresponds to the number of blocks found at beat S3, S4 or S5 for a 
given (S2, S3) pair. Sloan and SSmh, are as defined in the methods section. Total 
counts = 67 205. 

S2 stimulus. The S2 stimulus was delivered just after the refractory period of the last SI beat 
so that / at the site of stimulation was short (see figure 1). The short 7 produced slow V, which 
increased / for cells distal to the stimulus site. This initial dispersion in / following S2 was 
magnified by steep D restitution at each subsequent beat, ultimately leading to conduction block 
several beats later. 

Second, the value of the V restitution curve at 7 in j n was an important parameter in determining 
the likelihood of conduction block, as illustrated by figure 4. Panel (a) shows V restitution curves 
with four different values for the W^) (90, 75, 50, and 25% of V^). W^) in the wild type 
model was 75% of V^. Panel (b) shows a histogram for V (/^J = 50% of and panel 
(c) shows the histogram for 7(7,^) = 90% of V m . For V(/ min ) =. 90% of V^, the model 
had roughly 30% more instances of conduction block than in the wild type model, whereas the 
model with the lower cut-off values produced dramatically fewer instances of conduction block. 
In tact, the model with V^) = 25% of produced no cases of type II block. This result 
did not depend on the steepness of the V restitution curve. For example, no cases of type II block 
were found for a model with 0 = -1.4 and 5 = 2.0, which also had V(/ m m) = 25% of V^, 
but had a much steeper slope. The lack of type II block when V^) = 25% of can be 
understood by noting that if the V restitution curve approaches zero at 7,,,^, then a propagating 
wave that approaches a region with a very small 7 can slow down to allow 7 to increase just in 
front of it, permitting continued conduction. 

3.3. Contribution of D restitution to conduction block 

As mentioned above, D restitution contributed to conduction block by magnifying heterogeneity 
in 7. If the D restitution function had weak dependence on 7, any initial dispersion of 7 due to V 





Figure 4. The role of V restitution in type H block, (a) V restitution for four 
cases: 0 = 17.408 (black; wild type), 0 = 7.704 (red; V (/„*,) = 50% of 
VbJ, P = 2 028 (green; V(/ mia ) = 25% of V^) and = 30.236 (blue; 
V(/„nn) = 90% of V,™). (b) Histogram for V(/ mill ) = 50% of V**. Total 
counts = 8879 (cf 67 205 for wild type), (c) Histogram for V^) = 90% of 
Vnax- Total counts = 86032. There were no counts for V (/„*,) = 25% of Vmax- 
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Figure 5. The role of D restitution in type II block, (a) Steady-state D restitution 
for four cases: x D = 28 (black; wild type; max. slope = 1.09), x D = 32 (red; 
max. slope = 0.97), x D = 40 (green; max. slope = 0.82) and x D = 20 (blue; 
max. slope = 1.45). (b) Histogram for x D = 32. Total counts = 39 688. (c) 
Histogram for x D = 40. Total counts = 21 185. (d) Histogram for x D = 20. 
Total counts = 364 301. 

restitution would be lessened on subsequent stimulations. Figure 5 shows the effect of altering 
D restitution slope by increasing and decreasing the parameter x D . Panel (a) shows steady state 
restitution curves for four different values of x D . The wild type model was x D = 28. As shown 
in panels (b) and (c), increasing x D led to a significantly smaller peak in the histogram, as well 
as fewer total instances of conduction block. Conversely, decreasing x D , which led to a steeper 
steady state restitution slope, had significantly more counts (panel (d)). However, it is important 
to note that even for the x D = 40 model, which had a maximum steady state slope of 0.82, 
instances of type II block were found. 

3.4. Contribution of M to conduction block 

Because increasing M in the coupled maps model has been shown to eliminate steady state 
alternans, it was hypothesized that increasing M also would lead to fewer instances of type II 
block Figure 6 illustrates the effect of changing M on the likelihood of conduction block in the 
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model by increasing and decreasing the parameter a. Panel (a) shows steady-state restitution 
for four different values of a. The wild type model was a = 0.2. In panel (b) a = 0.58, which 
eliminated alternans but maintained a slope > 1. In fact, the steady state restitution slope was 
1 .62, even larger than in the wild type model. In panel (c) a = 0.8 and the maximum steady state 
restitution slope was 2.65. Both panels (b) and (c) show smaller peaks and fewer total counts 
than the wild type model. In panel (d) M was decreased (a = 0. 1) and the maximum steady state 
slope is 1 .09. The decreased memory model shows more instances of conduction block than the 
wild type model. This result did not depend on the magnitude of D. For example, shifting the 
curve from panel (c) by adding 70 ms to the function / (M B+ i , /„) did not produce more instances 
Of type II block (in fact slightly fewer (39 229) were found). Similarly, shifting the curve from 
panel (d) down by subtracting 70 ms from f(M^ u /„) produced roughly the same number of 
blocks (94 151) as in panel (d). Two other observations can be made. First, increasing M did 
not seem to be as effective as increasing z D in reducing the likelihood of block. Second, the 
steady-state D restitution function is not predictive of the likelihood of block. 

3.5. Contribution of electrotonic interactions to conduction block 

Electro tonic effects have been shown to prevent alternans in some models [7]. Therefore, 
changes in electrotonic current may be expected to have an effect on the incidence of conduction 
block. The coupled maps model allowed for the complete removal of electrotonic interactions by 
elirninating the diffusion terms in equation (4). As shown in figure 2, the potential for conduction 
block was present when the ^previous beat generated a short-to-long patternin/. This pattern led 
to a short-to-long pattern in the following D. In the wild type model, both the first and second 
derivative diffusion terms tended to decrease the long D and increase the short D, reducing 
the amount of heterogeneity and thereby reducing the likelihood of type II block. Accordingly, 
elimination of the diffusion terms was expected to increase the incidence of type II block and 
this expectation was confirmed; total counts of type H block were 88 603 (histogram not shown). 

4. Discussion 

In this study type II conduction block following the delivery of multiple premature stimuli was 
caused by the same sequence of events that caused conduction block during sustained rapid 
pacing in our previous study [12]. Stimulation at a short / produced a gradient of increasing 
/ (and, consequently, of D and V) along the fibre. The succeeding stimulus, if delivered at an 
appropriate interval, encountered a progressively decreasing / as it propagated down the fibre 
(i.e., the new wavefront encountered the waveback of the previous excitation). Thus, a short- 
long sequence of D at the site of stimulation was associated with the development of a long-short 
sequence of D at the opposite end of the fibre, similar to the pattern established during stable 
discordant alternans [19, 20, 24]. When the gradient in / along the fibre was sufficiently steep in 
the ascending direction, the subsequent action potential encountered an / < /„,&, and conduction 
block occurred. 

Previous studies have shown that the location of the block can be predicted, given the 
spatial profile of D for the preceding beat [22]. If this information is not known, however, then 
developing a quantitative theory for the type II block observed in this study becomes difficult, 
not only because the model consists of complicated nonlinear difference equations in both space 
and time, but also because the phenomenon of interest is a transient Therefore, the usual 
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Figure 6. The role of M in type II block, (a) Steady-state D restitution for 
four cases: or . != 0.2 (black; wild type; max. slope — 1.09), o = 0.58 (red; 
max. slope = 1.62), a = 0.8 (green; max. slope = 2.65) and or = 0.1 (blue; 
max. slope = 1.09). (b) Histogram for a = 0.58. Total counts = 48 387. (c) 
Histogram fora = 0.8. Total counts = 43715. (d) Histogram for or = 0.1. Total 
counts = 87433. 

tools of nonlinear dynamics, such as finding steady-state solutions and analysing their stability, 
do not apply. Still, inspection of the defining equations produces a qualitative understanding 
of the results presented in this study. The mechanism for conduction block relies on three 
determinants. First, a non-constant V restitution function produces spatial heterogeneity in 
7. Next, this heterogeneity is amplified (or at least it does not decay rapidly) due to sensitive 
dependence of Don/. Finally, the V restitution function does not drop sufficiently close to zero 
at very short /, causing the wave to collide with the refractory period of the preceding excitation. 

Preventing the first aspect of type II block could only be done by forcing the V restitution 
curve to be constant Even in that case, type II block could still occur if there were intrinsic 
heterogeneity in the fibre. Therefore, we explored how the likelihood of type II block could be 
diminished by two different methods. The first method was to alter the V restitution function so 
that conduction slowed at very short /, thus preventing encroachment on the refractory period. It 
is important to note that this method is only viable if the decrease in conduction velocity occurs at 



short/. Decreasing maximum conduction velocity isin fact pro-arrhythmic, became it decreases 
the wavelength [9]. Although selective slowing of conduction at short rest intervals may create a 
greater dispersion of refractoriness after a premature beat, the benefits of this intervention seem 
to outweigh the risks, at least under the conditions of this study. Whether a similar strategy will 
suppress conduction block in 2D or 3D tissue with realistic intrinsic heterogeneity remains to be 
determined. ... 

The second method for reducing the likelihood of type H block was to reduce the sensitivity 
of D on /. If the system is uni-dimensional, this is done by decreasing the D restitution slope. 
The lower the slope of the D restitution relation, the lower the likelihood of spatial dispersion 
of repolarization secondary to dispersion of /. However, the model used in this study did not 
have a simple relationship between the slope of the steady-state D restitution function and the 
likelihood of block. Instead, the important quantity is the sensitivity of D to the preceding I, 
given by 

v _!£ = M + iL!i (6) 

Y ~ dl 31 dMdl' 
The quantity y is the amplification factor relating the dispersion in D to the dispersion in the 
preceding /. If y is small, then any initial dispersion in I will decay quickly. If y is large, the 
dispersion of / can persist or even be amplified. It can be shown that at short cycle lengths the 
first term dominates, giving 

M* and /* are the steady state values at a given cycle length. This equation is of limited 
quantitative use since the system never reaches steady state (and Ihe cycle length is continuously 
changing). Still, it suggests that the likelihood of block can be reduced by either increasing x D or 
by increasing a. The simulation results support this qualitative argument This result highlights 
the fact that the dynamics during premature stimulation cannot be predicted from the steady-state 
D restitution function. In fact, the more relevant experimental measure would be the S1-S2 D 
restitution function at very short S2. The slope of this curve at short / would directly measure 
the dispersion of D after a very short S1-S2 interval had produced some initial dispersion in 
7. However, relating the S 1-S2 curve to any quantitative prediction about type H block is not 
straightforward. . 

Programmed electrical stimulation has been used extensively in the past to induce ventricular 
arrhythmias in patients thought to be at risk for the development of ventricular fibrillation [1]. 
However, the stimulation protocol used in the clinical studies differs significantly from that used 
here, in that arrhythmia induction was attempted using tightly coupled premature stimuli, which 
produced block similar to the type I block described in our study. Once the patient's anhythmia 
had been induced using multiple premature stimuli, the test was repeated after administration of 
an antiarrhythmic drag. In many cases a drug could be found that would suppress the patient's 
inducible arrhythmia. Yet the results of large scale clinical trials eventually revealed that patients 
sent home on such medications were as or more likely to die suddenly than patients who did 
not receive such treatment [3, 23]. Clearly, the induction of arrhythmias using this particular 
stimulation protocol was not predictive of the development of ventricular fibrillation. 

The failure of the standard method of programmed stimulation to accurately assess 
vulnerability to ventricular fibrillation may relate to the fact that this method was designed 
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primarily to induce conduction block and re-entry by perturbing intrinsic heterogeneity. The 
potential contribution of dynamic heterogeneity to the development of conduction block would 
not be assessed adequately using such a protocol. In addition, interventions that were intended to 
prevent the induction of arrhythmias using the standard protocol, such as slowing of conduction 
or prolongation of refractoriness, might be effective in that regard, secondary to a reduction in 
intrinsic heterogeneity, but might not reduce dynamic heterogeneity. Consequently, me beneficial 
effects of such interventions on arrhythmias induced by the standard stimulation protocol might 
not correspond to a reduction in the incidence of ventricular fibrillation, if induction of the latter 
is influenced importantly by dynamic heterogeneity. 

Dynamically induced heterogeneity and conduction block is one of many potential 
mechamsmsforspiraIbreak-upmcarmactissue(seeforexample[9]). Because wave propagation 
in the heart is extremely complex and influenced by many factors, it is unlikely that only 
one mechanism is responsible for induction and maintenance of all cardiac arrhythmias. 
Nevertheless, the results of this study suggest several potential interventions that may reduce 
the probability of aimythmiia induction, mcluding flattening the D restitution slope, decreasing 
V at short / and increasing M. There is experimental support for the idea that flattening D 
restitution suppresses ventricular fibrillation [21], but altering V and M with the objective of 
preventing ventricular fibrillation has not been attempted. Tests of these options 'may not be . 
likely in the near future, however, given that drugs that increase the recovery of sodium channels 
at depolarized potentials are not, to our knowledge, currently available, and alteration of M may 
not be possible until the ionic basis for this phenomenon has been established. Nevertheless, 
judicious alterations of D, V and M may hold future promise for new and more effective therapies 
for ventricular fibrillation. 
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Stimulation protocol for initiation of ventricular fibrillation; 3 methods for decreasing the risk of VP initiation. 



Ventricular fibrillation (VF), an extremely rapid and irregular heartbeat, is a leading cause of death in the Western world. Pharmacol gical 
approaches for preventing VF, which have included drugs that block sodium, calcium, and potassium channels, have been largely 
unsuccessful. Alternatively, implantable defibrillators can rescue a patient's heart from VF, but these devices are expensive and typically 
make use of a high-energy shock that can be painful and damaging to tissue. Thus, both methods for treating VF need improvement The 
development and testing of both treatments require a method for assessing the ability of a treatment to prevent or halt VF. Currently, VF can 
be initiated in the heart by a variety of methods, but none of these methods is similar to the clinical induction of VF, which typically occurs 
after a few (3- 10) premature ventricular activations. Recently, we developed a method for initiating VF using a computer model of normal 
one dimensional cardiac fibers. This model is based on the recovery properties of cardiac action potential duration and of the conduction 
velocity of electrical waves in the heart. The stimulus protocol consists of interrupting the normal sinus rhythm of the heart with 4 premature 
stimuli. If these stimuli are delivered with the appropriate timing (as predict by the model), the excitatory wave initiated by the fourth 
stimulus stops abruptly before reaching the end of the fiber. In two and three dimensional tissue, this "conduction block" can lead to wave 
break and the initiation of spiral wave reentry. Further, the model has shown that the likelihood of this conduction block occurring after the 
initiati n protocol is applied can be decreased by altering the recovery properties in three different ways: by altering the velocity recovery 
function so that conduction slows at short rest intervals, by altering the action potential duration recovery function to reduce the sensitivity 
of action potential duration to the preceding rest interval, and by increasing cardiac memory. Thus, this work has led to tw significant 
developments: 1) a stimulus protocol for initiating fibrillation in heart tissue that is similar to the clinical induction of VF and 2) three 
methods for altering the recovery properties of heart tissue to reduce the likelihood of VF induction. 
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A. Specific Aims 

The studies conducted thus far for this project 1 " 9 , coupled with the studies of several other 
groups 10 " 27 , have provided substantial evidence for a key role of electrical restitution in 
determining the dynamical behavior of cardiac tissue. These observations may be particularly 
relevant for understanding (and preventing) the development of ventricular fibrillation (VF). We 
now plan to extend our studies to more critically evaluate the contribution of dynamically- 
induced heterogeneity of repolarization in the setting of intrinsic electrical heterogeneity. To 
achieve these goals we propose the following specific aims: 

1. To determine the extent to which dynamic heterogeneity of repolarization 
predisposes to the development of conduction block and reentrant excitation. 

Our previous studies have shown that dynamic heterogeneity of repolarization, as induced 
by rapid pacing, predisposes to the development of conduction block in 1 -dimensional 
computer models and in isolated Purkinje fibers. We now propose to determine whether the 
delivery of 1-4 premature stimuli increases dynamic heterogeneity of repolarization and thereby 
precipitates conduction block. Studies will be conducted initially in computer models of a 
homogeneous 1 -dimensional cable and in isolated Purkinje fibers, where increased spatial 
dispersion of repolarization during premature stimulation can be assigned principally to the 
development of dynamical heterogeneity. 

2. To determine the extent to which dynamic heterogeneity of repolarization interacts 
with intrinsic heterogeneity in the development of conduction block and reentrant 
excitation. 

Studies will be conducted in experimental preparations and computer models having 
intrinsic spatial gradients of repolarization. Experimental models include arterially perfused 
canine left ventricle and anesthetized closed-chest German shepherd dogs having an inherited 
predisposition to ventricular arrhythmias and sudden cardiac death. 

The results of these studies will substantially improve our understanding of the 
mechanisms by which VF is initiated and sustained. As such, they promise to provide a basis 
for new and more effective therapies for the prevention of sudden cardiac death. 

B. Background and Significance 

Electrical restitution and spatial dispersion of repolarization 

Although several theories regarding the cellular electrophysiologic mechanism for VF 
currently are being evaluated 17,26,2 , substantial evidence has accumulated over the last few 
years to suggest that fibrillation is a state of spatiotemporal chaos consisting of the perpetual 
nucleation and disintegration of spiral waves 32,33 , in association with a period doubling 
bifurcation of local electrical properties 3,12,18,34 . Nucleation of the initiating spiral wave pair is 
caused by local conduction block (wave break) secondary to spatial heterogeneity of 
refractoriness in the ventricle 32,33,35 " 37 . Until recently, spatial heterogeneity was thought to result 
solely from regional variations of intrinsic cellular electrical properties 37, or from stimulation at 
more than one spatial location 21,25,36 . However, it is now appreciated that purely dynamical 
heterogeneity can be sufficient to cause conduction block during single site stimulation in both 
homogeneous one-dimensional models of canine heart tissue and in rapidly paced canine 
Purkinje fibers. A similar mechanism has been shown to precipitate conduction block and spiral 
break up in models of homogeneous two-dimensional tissue 11 . 
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The period doubling bifurcation implicated in the transition to conduction block is manifest 
as alternans, a beat-to-beat long-short alternation in the duration of the cardiac action potential 
3,12,18,25,34 p rev j OUS investigators have hypothesized that alternans can be accounted for by a 
simple uni-dimensional return map called the action potential duration restitution function 
10,14,39,40 Tnjs hypothesis assumes the duration of an action potential {APD) depends only on 
its preceding diastolic interval (D/) through some function f(DI) that is measured experimentally. 
If the APD restitution function has a slope > 1 , then a period doubling bifurcation occurs for 
some value of the pacing cycle (BCL), where BCL - APD + Dl. The velocity at which an action 
potential propagates (CV) also can be described by a restitution function, where CV = c(DI). 
More recently, other potential modulators of restitution, such as memory and calcium cycling, 
have been incorporated into formulations for restitution 6,11 41 . 

It has been shown previously that the combination of a steeply sloped APD restitution 
function and a monotonically increasing CV restitution function is sufficient to produce 
dynamical conduction block during sustained pacing at a short cycle length 1,42 . This 
observation may provide a generic mechanism for wave break and the onset of ventricular 
tachycardia and fibrillation. However, it is unlikely that the conditions used to demonstrate this 
phenomenon experimentally apply to the clinical situation, where the induction of ventricular 
tachyarrhythmias typically is associated with the interruption of normal cardiac rhythm by only a 
few premature beats. A single premature beat is sufficient to cause spatial heterogeneity in the 
form of discordant alternans 25 , but the conditions required for the development of conduction 
block in this setting have not been studied extensively. 

Therefore, we determined recently whether dynamic heterogeneity and conduction block 
occur in a computer model of canine heart fibers in which pacing at a slow rate is interrupted by 
1-4 premature stimuli 8 . This protocol simulates the interruption of sinus rhythm by 1-4 
premature ventricular complexes, a situation that can lead to the onset of ventricular fibrillation 
clinically. Furthermore, we determined whether the conduction block induced in this way can 
be accounted for by the same dynamical mechanism that underlies the development of 
conduction block during rapid pacing. A coupled maps model of heart tissue was used for the 
study, which allowed us to identify the roles that CV and APD restitution play in the 
development of dynamic heterogeneity. We also considered the additional contribution of 
cardiac memory (M) to the development of dynamical heterogeneity and conduction block. 

As described below, these studies indicated that a short-long-short-intermediate coupling 
interval pattern of premature stimuli induced marked spatial dispersion of repolarization and 
conduction block. The dynamical mechanism for the development of block was the same as for 
the development of discordant alternans during sustained rapid pacing. This type of conduction 
block could lead to wave break and spiral wave initiation if it occurred in two- and three- 
dimensional tissue, provided the block was local. The development of local block would be 
facilitated in intact myocardium by twist anisotropy and intrinsic heterogeneity. We have 
elected to pursue this idea, which provides the basis for the first specific aim of this project. 

An attractive experimental model for testing the hypothesis that maximizing dynamically- 
induced spatial heterogeneity of repolarization is conducive to the development of VF is the 
inherited sudden death model developed by Moi'se and co-workers 43 " 50 . This model consists of 
a colony of German shepherd dogs that are afflicted with an inherited propensity for sudden 
death. Affected dogs have frequent ventricular arrhythmias (VA), with those exhibiting 
ventricular tachycardia (VT) having the highest incidence of sudden death 43 . Death is 
unexpected and most frequent during sleep or at rest after exercise 44 . A window of vulnerability 
for the presence of VA and sudden death exists between 12 and 50 weeks of age with the peak 
affectedness between 20 and 28 weeks of age 46 . As determined by ambulatory Holter analysis, 
affected dogs can have up to 50,000 ventricular premature beats and runs of rapid VT during a 
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24 hr recording period. Two types of VT have been identified; (1) a rapid (heart rate > 300 
bpm), nonsustained, polymorphic form, and (2) a less rapid (heart rates < 250 bpm), sustained, 
monomorphic form 43,51, . The first type (85% of VT identified) is associated with sudden death 
and is frequently preceded by a pause. Moreover, this type of VT is most frequent during slow 
wave sleep, rapid eye movement (REM) sleep, and between 4 and 7 am 44 . The average 
incidence of sudden death in affected dogs is approximately 33%. However, the majority of 
affected dogs do not develop VF. The reason why most of these dogs do not develop VF, 
despite the presence of many episodes of rapid, polymorphic VT is unknown. In this regard, 
affected dogs are similar to some post-myocardial infarction patients, in that these patients may 
survive for years with a high incidence of VT before succumbing to sudden death 53 . 

The mechanism for the initiation and maintenance of VT in affected dogs has been 
investigated at several levels. The circumstances of death suggest that the autonomic nervous 
system is involved, in that the frequency of VT is greatest during rapid eye movement (REM) 
sleep 44 . During REM sleep 43 parasympathetic tone is high, but bursts of sympathetic tone also 
occur. Studies that examined sympathetic innervation to the heart using 123 l- 
metaiodobenzylguanidine scintigraphy and immunocytochemical staining for tyrosine 
hydroxylase revealed that large areas of the heart lack sympathetic innervation 54 . Further 
studies documented increased beta-receptor density, increased beta-adrenergic stimulation of 
adenylyl cyclase and abnormal beta-adrenergic signal transduction in the denervated regions of 
the heart 5 . In these same regions, cardiomyocytes are prone to the development of calcium 
oscillations, severe calcium overload, and cell death after exposure to f5-i-receptor agonists 55 . 

In addition to abnormalities of calcium cycling, other important alterations in the ionic 
currents responsible for repolarization are present in these dogs 52,56,57 , including a documented 
decrease in the density of the transient outward potassium current (l t0 ) and circumstantial 
evidence for decreased delayed rectifier (l Ks )- A reduction in outward repolarizing current 
presumably is responsible for the prolongation of action potential duration and development of 
prominent early afterdepolarizations (EADs) in left ventricular Purkinje fibers of affected dogs 50 . 
The EADs are the likely cause for the initiation of the pause dependent rapid polymorphic VT, 
the suspected initiator of VF. Abnormal T waves are seen on the surface ECG and are a 
reflection of the afterdepolarizations or the spatial inhomogeneity of repolarization among 
different regions of the left ventricle of affected dogs 49 . The latter could provide a permissive 
substrate for reentry arrhythmias as well. 

Thus, the German shepherd model provides an opportunity to determine whether certain 
patterns of multiple premature beats (secondary to EAD-induced triggered activity) are 
sufficient to induce VF in the setting of significant intrinsic heterogeneity of repolarization 
(secondary to repolarization alterations associated with sympathetic denervation). As 
discussed below under Preliminary Studies, we have begun to investigate whether patterns of 
premature stimuli designed to maximize spatial dispersion of repolarization via augmentation of 
dynamic heterogeneity initiate VF in these animals. Moreover, we also are testing whether the 
spontaneously occurring VT cycle lengths in these animals typically do not induce VF because 
they do not augment dynamically-induced spatial dispersion of repolarization sufficiently. 

C. Preliminary Studies 

Dynamically-induced spatial dispersion of repolarization 

Our previous theoretical and experimental studies of the impact of dynamic heterogeneity 
of repolarization on the development of conduction block were conducted in computer models 
of a homogenous 1 -dimensional cable and in isolated Purkinje fibers during fixed pacing. As 
the pacing cycle length was shortened progressively, a transition from 1:1 phase locking to a 
concordant alternans (2:2 locking) to discordant alternans to intermittent 2:1 conduction block 
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occurred 1 . These results confirmed the predictions of previous computer modeling studies of 
these phenomena 25 and extended experimental observations that the development of 
discordant alternans was associated with local conduction block and the initiation of reentry 19,59 . 
In addition, they demonstrated that complex dynamics was possible in a homogenous 1- 
dimensional cable, in the absence of intrinsic heterogeneity or anisotropy, and that the 
mechanism for the dynamics could be described using simple unidimensional maps for the 
restitution of APD and CV. 

The evolution of dynamics during fixed pacing occurs relatively gradually, which allows 
this phenomenon to be studied in the quasi-steady-state. As a result, a detailed description of 
the dynamics and identification of the underlying mechanism is possible. It is unlikely, 
however, that prolonged activation at a short constant cycle length initiates VF in patients at 
risk for sudden death. Therefore, we recently have begun to investigate a scenario that is more 
common in such patients - initiation of VF following a run of non-sustained VT. For these 
studies, computer models of a homogenous 1 -dimensional cable were generated using a 
coupled maps memory model 6 and an ionic model 2 . The models were paced at a cycle length 
of 1000 or 500 ms for 50 beats, after which a series of premature stimuli was delivered (Figures 
1 and 2). Premature stimuli delivered at short Dl (e.g., S 2 in Figure 2) elicited action potentials 
at the site of stimulation that, in accordance with the relevant APD and CV restitution relations, 
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Figure 1. Schematic drawings of: a. the 
coupled maps memory model, where action 
potential duration (D) is a function of the 
diastolic interval (I) and memory (M) during 
pacing at a fixed interval (T) and; b. the 
stimulation protocol to maximize spatial 
dispersion of repolarization. 




Figure 2. Spatial dispersion of repolarization induced in 
an ionic model of a homogenous 1 -dimensional cable by 
multiple premature stimuli (as in figure 3b). Note that 
short duration action potentials at the site of stimulation 
increase in duration as they propagate along the fiber, 
whereas long duration action potentials shorten. The 
resulting discordant alternans pattern culminates in 
conduction block. 



were short in duration and conducted slowly. Slow conduction along the fiber provided 
progressively more recovery time as the action potential conducted down the fiber. 
Accordingly, Dl increased along the fiber, as did the corresponding APD, resulting in a 
discordance between APD at the site of stimulation and at the distal end of the fiber. 
Conversely, premature stimuli delivered at longer Dl (S3 in Figure 2) produced longer duration 
action potentials that conducted more rapidly. The more rapidly conducting action potentials 
impinged progressively on the preceding Dl as they conducted along the fiber, creating a 
shorter Dl and corresponding APD. Once again, discordant APD between the site of 
stimulation and the distal end of the fiber resulted. The spatial dispersion of repolarization 
caused by this sequence of cycle lengths could be amplified by the delivery of multiple 
premature stimuli, eventually resulting in conduction block (S 5 in Figure 2). 
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To more completely characterize this sequence of events, we determined all possible 
combinations of S2-S5 that produced conduction block in the coupled maps model (block also 
occurred with many fewer combinations of S2-S3 and many more combinations of S 2 -S 7 ). 
Although approximately 67,000 combinations produced conduction block, the distribution of 
such combinations was not random, but was clustered according to specific S 2 and S 3 intervals, 
as shown in Figure 3. The incidence of block was greatest when S 2 was delivered at a coupling 
interval near SiS 2m in (the shortest S1S2 interval that produced a propagated response) and S 3 
was delivered at a coupling interval approximately 50 ms longer than S 2 S 3min . Although not 
shown on this 3-dimensional plot, the S 3 S 4 interval most often associated with block was near 
S 3 S 4m in, whereas the S 4 S 5 interval was approximately 10 ms longer than S 4 S5 mi n. This short- 
long-short-intermediate pattern of stimulation produced discordant alternans and conduction 

block of the type shown in Figure 2. In this 
study we also determined the impact of 
changes in the slope of the APD restitution 
relation, the magnitude of memory, the 
shape of the CV restitution relation and the 
strength of cell coupling. In general, the 
incidence of conduction block was 
reduced by decreasing the slope of the 
APD restitution relation, increasing 
memory, increasing cell coupling and 
extending the CV restitution relation to 
permit slow conduction (thereby 
preventing "head-tail" interactions) 8 . 

More recently, we have begun to 
investigate whether dynamically-induced 
spatial dispersion of repolarization 
promotes the development of conduction 
block and reentrant excitation in settings 
where underlying intrinsic heterogeneity is present. To that end, the sequence of premature 
stimulus intervals predicted by the computer model to cause conduction block was delivered to 
the right ventricles of closed chest anesthetized normal beagles and MAP recordings were 
obtained from the right and left ventricles. The dogs initially were paced at a cycle length of 
400 ms and a single premature stimulus (S 2 ) was delivered at progressively shorter SiS 2 
intervals until S 2 failed to capture (SiS 2mi n). The S-|S 2 interval was then set to SiS 2min + 40 ms 
and an S 3 stimulus was delivered at progressively shorter S 2 S 3 intervals until S 3min was 
encountered. The procedure was then repeated to find S 4min and S 5m in- Thereafter, a 
sequence of premature stimuli consisting of SiS 2mi n + 5 ms, S 2 S 3m in + 50 ms, S 3 S 4min + 5 ms 
and S 4 S 5min + 5-10 ms was delivered. As shown in Figure 4, the patterns of APD produced by 
this sequence of premature cycle lengths (CL V f) was similar to that produced in the computer 
models (c.f. Figure 2). Short duration responses (S) at the RV recording site were associated 
with longer duration responses (L) at the LV site and vice versa and the pattern of discordant 
alternans preceded the development of VF. 























Figure 3. Incidence of conduction block (counts - right 
axis) following the delivery of 4 premature stimuli (see 
protocol, lower left) in the coupled maps memory model. 
Results are shown for combinations of the S 1 S 2 and S 2 S 3 
intervals, referenced to S 2m i n and S 3mi n, respectively (S4 
and not shown). 
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Figure 4. Examples of MAP recordings obtained from the right 
(RV) and left (LV) ventricles of a closed-chest anesthetized 
beagle during delivery of multiple premature stimuli (CL VF ) to the 
RV. CL VF were chosen to produce maximum dynamically- 
induced dispersion of repolarization, as predicted from the results 
of the coupled maps model shown in Figure 3. 



The delivery of CL VF has 
produced VF in all 8 dogs tested 
thus far, whereas the delivery of 
multiple premature stimuli at 
sequentially shorter cycle lengths 
has failed to induce VF in any of the 
dogs (Figure 5). Although we have 
not, for obvious reasons, tested all 
67,000 possible combinations of 
cycle lengths predicted to cause 
conduction block in the model, nor 
have we tested more than a few of 
the intervals that are not expected 
to cause block, the highly reproducible induction of VF by CL VF lying near the peak of the plot 
shown in Figure 3 suggests that stimulus protocols designed to maximize dynamically-induced 

spatial heterogeneity of repolarization may 
reliably promote the induction of VF. The 
exact nature of the spatial dispersion 
produced by such a protocol will require 
more extensive mapping than can be 
accomplished by two MAP recordings, but 
the observation that CUf produces 
discordant alternans between at least 2 RV 
and LV recording sites suggests that the 
pattern of APD generated by the computer 
models may pertain to the intact heart. 

As a further test of the hypothesis that 
dynamically-induced spatial heterogeneity is 
arrhythmogenic, we determined whether a 
reduction in the slope of the APD restitution 
relation suppresses the development of 
spatial heterogeneity APD and conduction 
block during premature stimulation, as 
shown in the computer model 8 . For these 
studies, we tested the effects of verapamil, which is known to decrease the slope of the 
dynamic restitution relation 5 , in the anesthetized closed chest dogs described above (n = 8). 
MAP duration (MAPD) was measured in RV and LV endocardium during programmed 
stimulation of the RV using 4 appropriately timed premature stimuli (S 2 -S 5 ), following a train of 
20 stimuli at a constant S1S1 interval of 400 ms. The stimulation protocol incorporated the 
same short-long-short-intermediate sequence shown in Figure 5, adjusted as necessary to 
account for changes in S 2 min-S 5 min produced by verapamil. VF induction was attempted during 
control and 30 minutes after i.v. verapamil (0.1, 0.3 or 1.0 mg/kg/min). 

VF induction occurred in all dogs during control and after 0.1 mg/kg verapamil and was 
associated with discordant alternans of MAPD between RV and LV. In contrast, after the two 
higher verapamil doses MAPD alternans was reduced in magnitude and was concordant and 
VF could not be induced. These effects were not associated with significant changes in the 
slope of the standard (SiS 2 ) MAPD restitution relation, which was 0.41 +0.17 during control 
and was 0.47 + 0.28, 0.56 + 0.20 and 0.39 + 0.18 after verapamil. Thus, suppression of spatial 
APD heterogeneity by verapamil was associated with an inability to induce VF using premature 




Figure 5. Incidence of conduction block following the 
delivery of 4 premature stimuli in the coupled maps 
memory model (Figure 3), and following the delivery of 
CL VF in anesthetized beagles (filled circles). CUf 
induced VF in all 8 dogs. In contrast, patterns of cycle 
lengths that fell outside CL VF did not induce VF 
(unfilled circles). 
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stimuli. However, the doses, of verapamil required to produce this effect had significant 
negative dromotropic, chronotropic and inotropic effects: verapamil increased the PR interval 
(from 84.7 + 3.8 during control to 111.8 + 12.2, 132.0 + 15.0 and 139.5 + 18.0 ms after 
verapamil), decreased heart rate (from 139.0 + 9.7 to 110.0 + 2.1, 101.6 + 3.8, and 95.2 + 8.8 
bpm) and decreased mean blood pressure (from 116.8 + 7.3 to 107.4 + 11.4, 108.0 + 6.3 and 
101.2 + 9.4 mmHg). These results suggest that suppression of VT-induced spatial electrical 
heterogeneity may be one means of preventing VF, but that verapamil is unlikely to be the drug 
of choice for producing such an effect in patients. 

Although the induction of VF by CL V f in normal beagles confirms the predictions made 
by the computer models, this result is somewhat expected, in that the data used to construct 
the models was obtained from normal canine myocardium. This model has the additional 
limitation of requiring external stimuli to induce VF, since these dogs do not have spontaneous 
ectopy. To test the predictions from the computer model more rigorously, we next investigated 
whether CL V f induced VF in German shepherd dogs that display non-sustained polymorphic VT 
and have an inherited predisposition to sudden death (Figure 6). As discussed above, one of 
the interesting and unanswered questions with respect to these animals is why they often live 
for months with multiple episodes of rapid, polymorphic VT before succumbing to sudden 

death. Based on the response of the 
beagles to CL V f. it seems possible that the 
German shepherds, despite having copious 
amounts of ectopy, do not die sooner 
because they do not routinely generate the 
proper sequence of coupling intervals, 
insofar as promoting dynamically-induced 
spatial dispersion of repolarization is 
concerned. 

As an initial test of this idea, 5 affected 
German shepherd dogs were anesthetized 
and instrumented as described above for 
the beagles. CL V f were determined from 
measurements of S2min"S5min and 
predictions from the beagle-based 
computer model and were delivered to the 
RV or LV. VF was induced in 3 of the dogs, 
but not in the other two, even after several adjustments of the CL V f- Sequences of cycle 
lengths obtained from Holter recordings during multiple episodes of polymorphic VT (CLvt) also 
were delivered to the dogs. None of these sequences induced VF. The sequences during 
CLvt consisted of progressively shorter cycle lengths, in contrast to CL V f, which were short- 
long-short-short. In addition, even during longer runs of non-sustained VT, such as the run 
shown in Figure 6, a CL V f - type pattern did not occur. This observation raises the question of 
why do the German shepherds eventually die suddenly? Unfortunately, despite the hundreds 
of Holter recordings obtained since this colony was established, VF has never been captured. 
To remedy this situation, we currently are implanting Reveal Plus© (Medtronic) event recorders 
in severely affected dogs. Thus far, we have captured the development of VF in one dog. The 
sequence of cycle lengths preceding VF was similar to that used to precipitate VF in the 
programmed stimulation studies, but clearly more episodes of spontaneous VF need to be 
captured and analyzed before definitive conclusions can be drawn. 

D. Research Design and Methods 
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Figure 6. Spontaneous non-sustained polymorphic 
ventricular tachycardia in a German shepherd dog. 
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Specific Aim #1. 
Experimental plan 

The mechanism for conduction block following the delivery of multiple premature stimuli will 
be studied initially using a one-dimensional computer model of a canine heart fiber, as 
described previously 1 . The fiber will be stimulated at one end at a cycle length near normal 
canine sinus rhythm (S-|Si=500 ms). After 10 Si stimuli, 4 premature stimuli (S2, S3, S 4 and S 5 ) 
will be delivered. The coupling intervals of the premature stimuli will be varied and conduction 
of the resultant responses observed. For each combination of premature stimuli, one of three 
outcomes is possible: 1 ) the stimulus elicits an action potential that propagates along the entire 
fiber; 2) the stimulus does not elicit a propagated response (Type I block); 3) the stimulus elicits 
an action potential that blocks before reaching the end of the fiber (Type II block). 

To assess the vulnerability of the simulated tissue to Type II block, the S2S3S4S5 
combinations will be applied as follows: the S-iS 2 interval will be varied from the minimum value 
that conducts (S 1 S 2m in) to SiS 2m in + 20 ms. For each S-|S 2 interval, the S 2 S 3 interval and the 
S3S4 interval will be varied in combination from the minimum value that conducts up to a value 
of 250 ms for each interval. For each S 2 S 3 S 4 combination that conducts, the S 4 S 5 interval will 
be varied from the minimum interval that generates an action potential at the site of stimulation 
to the minimum interval that generates an action potential that conducts down the entire fiber. 
If no S 2 S 3 S 4 S 5 combination is found that produces conduction block for more than a 20 ms 
window in either the S 2 S 3 interval or the S 3 S 4 interval, the search is halted for that interval. 
Intervals are varied in steps of 1 ms. 

The computer model currently is based on the coupled maps memory model, which may 
not replicate certain aspects of the dynamics during premature stimulation 41 . Accordingly, 
future versions of the model will incorporate the results obtained from the standard restitution 
protocol for each of the four premature stimuli. The results generated by the latter model will be 
compared with the results generated by the current model and will be evaluated with respect to 
the success or failure of the CL V f predicted by the different models to induce VF in the arterially 
perfused left ventricular preparations and in the whole animal experiments (see below). Other 
models of restitution 11,41 also may be tested, as warranted. The restitution relations for APD 
and CV in the current model were derived from data obtained from canine endocardium. In 
future versions of the model, these functions will be obtained from epicardial and mid- 
myocardial layers of myocardium as well. 

Once the effects of multiple premature stimuli on spatial dispersion of repolarization have 
been determined in the homogenous computer model, simulations will be conducted in the 
coupled maps model after introduction of a spatial gradient of intrinsic repolarization, using a 
modified version of the method of Henry and Rappel . The repolarization gradient and pattern 
of cell coupling will be designed to mimic the transmural gradient of repolarization reported for 
canine myocardium 61 " 63 . The interactions between intrinsic and dynamic heterogeneity of 
repolarization also will be studied in 2- and 3- dimensional computer models, as they become 
available as the result of studies to be conducted in collaboration with Drs. Jeffrey Fox and 
Wouter-Jan Rappel ("Computer model of the canine ventricle", supported by HL075515). 

Once the results of the computer simulations are known, experiments will be conducted in 
isolated canine Purkinje fibers. These experiments should, based on our previous experience 1 , 
be an appropriate test of the effects of multiple premature stimuli on spatial dispersion of 
repolarization in a (nearly) homogenous, (functionally) 1 -dimensional cable. Thereafter, the 
effects of CLvf will be tested in isolated arterially perfused canine left ventricle, a preparation 
capable of generating VF 5 . Tests of the CL V f protocol also will be conducted in normal beagle 
dogs, where the results generated by the computer model, which is based on data obtained 



R746281.I 



30 



from such hearts, is expected to be relevant (Figures 4 and 5). Finally, the CL V f protocol will be 
tested in German shepherd dogs having an inherited predisposition to sudden death. 
The details of the various experimental protocols are given in the following sections. 

One-dimensional computer models 

Our initial studies will be conducted using a coupled maps model of a one-dimensional 
cardiac fiber, as described previously 1,94 . The model is based on the equation 

r.«(x,)=T.A*,)-D IHrl (x i ). (1) 

T n+ i (xj is the time interval between activations of site x,. It is determined by including the time 
delays caused by the propagation from the pacing site to site x, , which yields 

x is the time interval between activations applied to the pacing site and Ax=0. 1 is the length of 
a single cell (time units in ms and space units in mm). The conduction velocity V n (Xi) depends 
only on / through the velocity recovery function V n = c(l n ) given by c(l)=V max (1 -exp(-(l+fi)/S)). 
V max =0. 72, S=14, and p is varied to adjust the value of V at l=l min . D is determined locally based 
on a memory model mapping 94 given by: 

M„ + , = g(M n , I„,D n ) = e'° ' r - [l + (M n -lK°" /r "] 

r B V ( 3 ) 

D n+l = /(M fl+1 ,/J = (1 - oM n+l )[A + l + eHln . C)/tD J 

x m is the time constant for memory. A = 88, B = 122, C = 40, and r m = 180 94 . a determines the 
influence of memory on D, and x D is used to adjust the dependence of the D recovery function 
on /. Coupling between sites is included by using the diffusion terms from Echebarria and 
Karma 95 . Including diffusion then yields 

D n+ ,=f(M ll+v I n ) + ?V 2 D n+l -wVD n+l , (4) 

with £ = 1.0 and w = .35. Discretizing the derivatives produces a tri-diagonal linear system of 
equations that can then be solved easily. The defining equation for the model is therefore 

M„ +1 =g(M tt ,I n ,D n ) 

D n+] =AM n+i ,I n ) + { 2 V 2 D n+ , -wVD„ +] (5) 

Cellular action potentials and wave propagation will be simulated using a modified version 
the Winslow ionic model, 66 as described in detail 2 ' 67 . Briefly, compared to the native Winslow 
model, l K i is decreased at depolarized potentials. 56 The maximum conductance and 
rectification of l Kr are increased and activation kinetics are slowed. l Ks is increased in 
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magnitude and activation shifted to less positive voltages. 68 L-type calcium current is modified 
to produce a smaller, rapidly inactivating current. Finally, a simplified form of intracellular 
calcium dynamics is used 69 The PDE for a one-dimensional piece of tissue will be solved 
numerically, as described by Qu et a/. 20 The length of the one-dimensional cable is set to 220 
cells (30 cells longer than the minimum length required for the development of at least one 
node), with the length of each cell = 200 ^M, and a diffusion coefficient of .001 cm 2 /msec. 

Isolated Purkinje fibers 

Adult beagle or mongrel dogs weighing 10-30 kg will be anesthetized with Fatal-Plus (86 
mg/kg i.v.) and their hearts will be excised rapidly 5 . Free running Purkinje fibers 1.5-2.5 cm in 
length and 1-2 mm in width are excised from either ventricle, mounted in a Plexiglas chamber 
and superfused with Tyrode solution at 15 ml/min. The Tyrode solution is bubbled with 95% 0 2 
and 5% C0 2 . The P0 2 is 400-600 mmHg, the pH is 7.35 + 0.05 and the temperature is 37.0 + 
0.5 °C. The composition of the Tyrode (in mM) is: MgCI 2 0.5, NaH 2 P0 4 0.9, CaCI 2 2.0, NaCI 
137.0, NaHC0 3 24.0, KCI 4.0 and glucose 5.5. The fibers are stimulated using rectangular 
pulses of 2 ms duration and 2-3 times the diastolic threshold (0.1-0.3 mA), delivered through 
Teflon-coated bipolar silver electrodes. Transmembrane action potentials recordings are 
obtained simultaneously from 4-6 sites along the fiber. The recordings are sampled at 2 kHz 
with 12-bit resolution using custom data acquisition programs written in Linux and run on a PC. 

The dependence of APD on the preceding Dl will be determined using a standard 
restitution protocol and a dynamic protocol, as described previously 9 . For the standard 
protocol, single test pulses (S 2 ) are delivered after every 20 th basic pulse (Si) at BCL = 400 ms. 
The SiS 2 coupling interval is progressively shortened in steps of 10-20 ms starting from 300 ms 
until the premature stimulus fails to capture. The duration of the action potential elicited by S 2 is 
measured at 95% of repolarization (APD95) and is plotted as a function of the preceding Dl, 
where the Dl equals the SiS 2 interval minus APD95 of the response to the last Si stimulus. This 
procedure will be repeated to measure restitution following delivery of an S 3 , S 4 and S 5 
stimulus, to generate data comparable to those generated in the whole dog experiments. 

The relationship between APD and Dl also will be determined during pacing at a constant 
BCL (dynamic restitution). The BCL is shortened from 400 to 200 ms in steps of 50 ms and 
from 200 ms to the effective refractory period in steps of 5-10 ms. At BCL that produce a 1:1 
stimulus.response locking, pacing is stopped after steady state has been reached and APD 95 of 
the last paced action potential is measured. During APD alternans, pacing is interrupted twice 
to directly measure APD 95 of both the long and the short action potential. 

The restitution data will be entered into the coupled maps computer model to generate 
sequences of multiple premature stimuli that maximize dynamically-induced spatial dispersion 
of repolarization and precipitate conduction block. The model-generated cycle lengths will then 
be delivered to the Purkinje fiber and the spatial patterns of APD and the incidence of 
conduction block determined. Cycle length sequences generated from both the standard and 
dynamic restitution relations will be tested to determine which sequence produces spatial 
patterns in the fiber that most closely agree with those predicted by the model. If the spatial 
patterns produced by both of these sequences differs significantly from that predicted by the 
model, the coupled maps model will be modified by altering the configuration for memory. 

Specific aim # 2 

Arteriallv perfused left ventricle 

Adult beagle or mongrel dogs or juvenile affected German shepherd dogs will be 
anesthetized as described above and their hearts will be excised rapidly and placed in cool 
Tyrode solution 5 . The circumflex coronary artery or a branch of the right coronary artery will be 
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cannulated using PE tubing and a transmural section of tissue measuring 30-50 mm in width, 
30-90 mm in length and 10-18 mm in depth excised. The preparation is suspended in a 
Plexiglas chamber with the epicardial or endocardial surface facing up or is placed on-end, to 
allow recordings to be obtained from both surfaces simultaneously. The preparation is 
perfused via the coronary artery and is superfused with normal Tyrode solution at a flow rate of 
35 ml/min. Perfusion pressure is 50-80 mm Hg and the temperature is 37.0 - 38.0°C. 

Electrical activity will be mapped using arrays of 64 monophasic action potential (MAP)-type 
recording electrodes, supplemented by 1-4 floating glass microelectrodes. In those 
experiments for which transmural recordings are required, 2-3 roving MAP electrodes or 
floating microelectrodes will be lowered onto the uppermost transmural edge of the preparation. 
The MAP and transmembrane action potential recordings are displayed on a storage 
oscilloscope and a computer monitor and are sampled at 1250 Hz with 12-bit resolution. The 
MAP signals are high-pass (cutoff = 0.15 Hz) and low-pass filtered (cutoff = 600 Hz). Records 
of 4-7 seconds duration are obtained every 20-40 seconds during the course of the experiment. 

The standard and dynamic restitution relations will be determined from microelectrode or 
MAP recordings obtained from epicardium, endocardium and midmyocardium, with the latter 
being accessed from the exposed transmural edge of the preparation. For epicardial and 
endocardial measurements, action potentials will be recorded from several sites from base to 
apex, with the expectation that a gradient of APD and restitution kinetics will exist. Following 
measurement of the restitution relations, multiple premature stimuli will be delivered to the 
preparations and their effect on spatial dispersion of repolarization, the development of 
conduction block and the induction of VF will be determined. 

Intact beagle dogs 

We plan to complete our studies using normal beagle dogs and have provided for such 
studies in our application for a Scientist Development Grant from the National American Heart 
Association. However, because the term of this application is 3 years, compared with a term of 
4 years for the national award, we have elected to not include the studies of beagle dogs in this 
application, so that Dr. Gelzer may move forward sooner with the other studies. If this award is 
funded, the beagle studies will be competed by a cardiology resident, with assistance and 
financial support from Drs. Moise and Gilmour, using other resources. 

Intact German shepherd dogs 

Studies will be conducted in German shepherd dogs having frequent ventricular premature 
complexes and non-sustained ventricular tachycardia (VT). Surface ECG recordings will be 
acquired using 24 hr ambulatory Holter monitoring. All data will be stored on 350 MB 
removable PC flashcards and automated Holter analysis will be performed using the Vision 
Premier™ software (Spacelabs Burdick Inc.). In each dog a set of 25 coupling intervals of 
spontaneous runs of VT (first 5 successive RR intervals of several episodes) will be measured 
to generate a library of coupling intervals (CLvt) to be used for programmed stimulation. On 
the day of the study, the dogs will be induced using pentobarbital (20 mg/kg IV) and maintained 
with a continuous rate infusion of pentobarbital (5 mg/kg/hr) and fentanyl dihydrogen citrate 
(0.004 mg/kg/hr). Introducer sheets will be placed percutanously into the right femoral vein and 
artery and MAP recording catheters will be advanced into the RV and LV apex. These 
catheters are used both for cardiac pacing and for MAP recording. The MAPs and a standard 
surface ECG will be recorded using a multi-channel physiologic recorder (Biopac Systems, Inc). 

The APD restitution relation will be determined using a dynamic restitution protocol, as 
described previously 9 . Briefly, the RV or LV will be paced at S-^ = 400 ms for 50 beats and 
then at progressively shorter S1S1 intervals. The APD restitution relation will be determined by 
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plotting APD at 90% of repolarization (APDg 0 ) of the last paced beat at each BCL as a function 
of the preceding Dl. APD restitution also will be determined using a standard protocol, where 
the RV or LV will be paced at S1S1 = 400 ms for 20 beats and then at progressively shorter 
S1S2 intervals, until the S2 stimulus fails to generate a propagated response. The S1S2 interval 
will then be set to 40 ms longer than the longest S1S2 interval that failed to capture and an S 3 
stimulus will be delivered at progressively shorter S2S3 intervals, until the S 3 stimulus fails to 
capture. This procedure will then be repeated using an S 4 and an S 5 stimulus. 

Once the APD restitution relations have been determined, the RV or LV will be paced at 
S1S1 = 400 ms for 20 beats, followed by a specific set (CI_vt) of premature stimuli (S2S3S4S5) 
obtained from the library of coupling intervals that occurred during spontaneous VT. The RV or 
LV will then be paced at S1S1 = 400 ms for 20 beats, followed by a specific set (CL V f) of 
premature stimuli (S2S3S4S5) that are expected, based on calculations using the computer 
model, to produce conduction block and VF. If VF occurs, the dogs will be converted to sinus 
rhythm using external cardioversion/defibrillation. 

Limitations 

A potential limitation of the studies proposed for this project is that the spatial complexity of 
repolarization in the intact heart will not be recreated in the computer models. However, 
previous studies in which a simple computer model has been used to predict the behavior of 
actual myocardium have been remarkably successful not only in predicting, but also in 
providing a mechanistic explanation for, complex dynamical behavior related to the initiation 
and perpetuation of ventricular tachyarrhythmias 1 . Consequently, the studies outlined in this 
proposal are likely to provide important, robust mechanistic insights that can be refined as 
necessary at a latter date when more realisitc computer models of the heart become available. 

A related limitation applies to the proposed studies of isolated cardiac tissue and intact 
hearts. In these studies, the complexity of actual cardiac tissue will be present, but the 
mapping techniques to be used will not be capable of fully evaluating that complexity. To do 
so would require simultaneous measurements of electrical activty from all sites, or at least from 
representative sites, throughout the heart. Nevertheless, even though the full range of spatial 
dynamics may not be captured in these studies, the proposed experiments are highly likely to 
indicate whether dynamically-induced heterogeneity of repolarization is an important 
determinant for the induction of VF. In the future, we anticipate that we will have the capability 
of recording electrical activity transmurally using nanofabricated multielectrode arrays. These 
devices currently are under development and support for this project has been requested. (NIH 
R01 HL073644-01, "MEMs sensors for arrhythmia detection and intervention", R.F. Gilmour 
Jr., PI, E. Bodenschatz, A. Lai, B. Lerman and D. C. Christini, Co-Investigators). 

Although we will evaluate the predictions from the computer models regarding which 
combinations of premature stimuli are likely to induce VF using isolated myocardium and intact 
hearts, it seems likely that the electrophysiological substrates in those preparations may differ 
significantly from those present in patients with ischemic heart disease or various 
cardiomyopathies. Consequently, the proposed studies will need to be extended in the future 
to experimental circumstances that mimic the substrates in patients at risk for sudden cardiac 
death. In that regard, we envision future studies in acutely ischemic myocardium, similar to 
those we have conducted previously, and in a pacing-induced heart failure model, where 
recent studies have indicated that augmented intrinsic dispersion of repolarization provides a 
permissive substrate for the development of reentrant arrhythmias 108 . 

A final limitation (or caution) relates to the overall hypothesis that steep restitution of APD 
predisposes to dynamically-induced dispersion of repolarization following premature 
stimulation and that the latter can be used to establish the risk for development of VF. 
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Programmed stimulation to identify patients at risk for sudden death and to evaluate the 
efficacy of therapy to prevent sudden death has, of course, been tried, with very disappointing 
outcomes 109 . Consequently, the results of the proposed studies need to be interpreted in light 
of the realization that cellular electrical dynamics is unlikely to be predicted with absolute 
fidelity by simple models and that provocative tests that yield useful information in one setting 
may not necessarily be useful in a different setting. In that regard, the mode of programmed 
stimulation used in CAST and other such studies was designed primarily to probe intrinsic 
heterogeneity of refractoriness. The studies outlined in this proposal are designed to provoke 
and to evaluate dynamical heterogeneity, alone or superimposed on intrinsic heterogenity. It 
seems likely that although our preliminary experiments in the normal beagles indicate that 
patterns of stimulation based solely on dynamical heterogenity reliably induce VF, intrinsic 
heterogeneity may need to be considered more fully in order for useful predictions to be 
generated for diseased myocardium. 

Ethical aspects of the proposed research 

Adult beagles or German shepherd dogs will be anesthetized with Fatal-Plus (390 mg/ml 
pentobarbital sodium; Vortech Pharmaceuticals; 86 mg/kg IV). Hearts will be removed via a 
thoracotomy and segments of the heart placed in a tissue bath for further study. All 
experiments have been approved by the Institutional Animal Care and Use Committee of the 
Center for Research Animal Resources at Cornell University. The animals will be housed at the 
AAALAC accredited Laboratory Animal Services facility at the College of Veterinary Medicine, 
Cornell University. Excellent support service and veterinary care are provided. 

The use of animal tissue is necessary because there is no currently available mathematical 
or computer model that adequately simulates the complex three-dimensional environment 
associated with the initiation and perpetuation of the type of heart rhythm disorders we plan to 
study. Canine hearts are used because their cardiac electrophysiological properties are similar 
to those of many other species, including humans, and their hearts are large enough to permit 
initiation and recording of rhythm disturbances. In addition, canine hearts have been studied 
extensively by other investigators, which facilitates data comparisons. 
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WHAT IS CLAIMED: 

1 . A method of evaluating the effect of a physiological condition on 
the occurrence of ventricular fibrillation, said method comprising: 

providing a test system; 

initiating a ventricular fibrillation inducing sequence in the test system by 
interrupting normal sinus heart rhythm with premature electrical stimuli; 

initiating a ventricular fibrillation recovery sequence in the test system 
following said initiating a ventricular fibrillation inducing sequence; 

subjecting the test system to a physiological condition to be tested before, 
during, or after said initiating a ventricular fibrillation recovery sequence; and 

identifying physiological conditions which affect ventricular fibrillation 
recovery in the test system. 

2. The method according to claim 1 , wherein the said physiological 
condition to be tested is velocity restitution. 

3. The method according to claim 1 , wherein the said physiological 
condition to be tested is action potential duration restitution. 

4. The method according to claim 1 , wherein the said physiological 
condition to be tested is cardiac memory. y 
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5. The method according to claim 1 , wherein said initiating a 
ventricular fibrillation inducing sequence is carried out by initiating 4 premature stimuli 
under conditions effective to initiate an excitatory wave. 

6. The method according to claim 1 , wherein said initiating a 
ventricular fibrillation recovery sequence is carried out by altering velocity recovery 
function to slow conduction at shorter rest intervals. 

7. The method according to claim 1 , wherein said initiating a 
ventricular fibrillation recovery sequence is carried out by altering action potential 
duration recovery function to reduce sensitivity of action potential duration to a preceding 
rest interval. 

8. The method according to claim 1 , wherein said initiating a 
ventricular fibrillation recovery sequence is carried out by increasing cardiac memory. 

9. The method according to claim 1 , wherein the test system is 
selected from the group consisting of a test animal, a tissue, a cell culture, or an in vitro 
system. 

1 0. The method according to claim 9, wherein the test system is a test 

animal. 
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11. A method of identifying treatment candidates as therapeutic 
strategies for prevention of ventricular fibrillation, said method comprising: 

providing a test system; 

initiating a ventricular fibrillation inducing sequence in the test system by 
interrupting normal sinus heart rhythm with premature electrical stimuli before, during, 
or after administering the treatment candidate to the test system; and 

identifying treatment candidates which prevent said initiating from 
inducing ventricular fibrillation as therapeutic strategies for prevention of ventricular 
fibrillation. 

1 2. The method according to claim 1 1 , wherein the treatment 
candidate is a pharmaceutical compound. 

13. The method according to claim 1 2, wherein the pharmaceutical 
compound is a calcium channel antagonist. 

14. The method according to claim 1 1 , wherein the treatment 
candidate is one or more electrical impulses. 

1 5. The method according to claim 1 1 , wherein said identifying 
identifies treatment candidates which achieve velocity restitution values consistent with 
histograms b and c of Figure 4 in "Dynamic Mechanism for Conduction Block in Heart 
Tissue" as therapeutic strategies for prevention of ventricular fibrillation. 
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1 6. The method according to claim 1 1 , wherein said identifying 
identifies treatment candidates which achieve potential duration restitution values 
consistent with histograms b and c of Figure 5 in "Dynamic Mechanism for Conduction 
Block in Heart Tissue" as therapeutic strategies for prevention of ventricular fibrillation. 

1 7. The method according to claim 1 1 , said identifying identifies 
treatment candidates which achieve cardiac memory values consistent with histograms b 
and c of Figure 6 in "Dynamic Mechanism for Conduction Block in Heart Tissue" as 
therapeutic strategies for prevention of ventricular fibrillation. 

1 8. The method according to claim 1 1 , wherein the test system is 
selected from the group consisting of a test animal, a tissue, a cell culture, or an in vitro 
system. 

1 9. The method according to claim 1 8, wherein the test system is a test 

animal. 

20. The method according to claim 1 1 , wherein said initiating a 
ventricular fibrillation inducing sequence is carried out by initiating 4 premature stimuli 
under conditions effective to initiate an excitatory wave. 

21. A method of identifying treatment candidates as therapeutic 
strategies for treating ventricular fibrillation, said method comprising: 
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providing a test system; 

initiating a ventricular fibrillation inducing sequence in the test system by 
interrupting normal sinus heart rhythm with premature electrical stimuli resulting in 
ventricular fibrillation in the test system; 

administering the treatment candidate to the test system undergoing 
ventricular fibrillation; and 

identifying treatment candidates which modulate ventricular fibrillation as 
therapeutic strategies for treatment of ventricular fibrillation. 

22. The method according to claim 2 1 , wherein the treatment 
candidate is a pharmaceutical compound. 

23. The method according to claim 22, wherein the pharmaceutical 
compound is a calcium channel antagonist. 

24. The method according to claim 21 , wherein the treatment 
candidate is one or more electrical impulses. 

25. The method according to claim 2 1 , wherein said identifying 
identifies treatment candidates which achieve velocity restitution values consistent with 
histograms b and c of Figure 4 in "Dynamic Mechanism for Conduction Block in Heart 
Tissue" as therapeutic strategies for prevention of ventricular fibrillation. 
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26. The method according to claim 2 1 , wherein said identifying 
identifies treatment candidates which achieve potential duration restitution values 
consistent with histograms b and c of Figure 5 in "Dynamic Mechanism for Conduction 
Block in Heart Tissue" as therapeutic strategies for prevention of ventricular fibrillation. 

27. The method according to claim 2 1 , said identifying identifies 
treatment candidates which achieve cardiac memory values consistent with histograms b 
and c of Figure 6 in "Dynamic Mechanism for Conduction Block in Heart Tissue" as 
therapeutic strategies for prevention of ventricular fibrillation. 

28. The method according to claim 21 , wherein the test system is 
selected from the group consisting of a test animal, a tissue, a cell culture, or an in vitro 
system. 

29. The method according to claim 28, wherein the test system is a test 

animal. 

30. The method according to claim 2 1 , wherein said initiating a 
ventricular fibrillation inducing sequence is carried out by initiating 4 premature stimuli 
under conditions effective to initiate an excitatory wave. 
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31. A method for evaluating the predisposition of a test animal for the 
induction of ventricular fibrillation from a condition of ventricular tachycardia, said 
method comprising: 

providing a test animal in ventricular tachycardia; 

monitoring the electrical stimuli in the heart of the test animal; 

determining if each stimulus in groups of 4 stimuli correspond to rest 
interval values predicted to lead to ventricular fibrillation consistent with the histogram of 
Figure 3 in "Dynamic Mechanism for Conduction Block in Heart Tissue." 

32. A method of identifying treatment candidates as therapeutic 
strategies for preventing ventricular tachycardia from developing into ventricular 
fibrillation, said method comprising: 

providing a test animal in ventricular tachycardia; 

monitoring the electrical stimuli in the heart of the test animal; 

determining if an initial 3 stimuli in groups of 4 stimuli in the heart 
correspond to rest interval values predicted to lead to ventricular fibrillation consistent 
with the histogram of Figure 3 in "Dynamic Mechanism for Conduction Block in Heart 
Tissue"; 

identifying treatment candidates which prevent occurrence of a fourth 
stimuli corresponding to a rest interval value predicted to lead to ventricular fibrillation 
consistent with the histogram of Figure 3 in "Dynamic Mechanism for Conduction Block 
in Heart Tissue"; 
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identifying treatment candidates which prevent said ventricular 
tachycardia from becoming ventricular fibrillation as therapeutic strategies for prevention 
of ventricular fibrillation. 

33. The method according to claim 32, wherein the treatment 
candidate is a pharmaceutical compound. 

34. The method according to claim 33, wherein the pharmaceutical 
compound is a calcium channel antagonist. 

35. The method according to claim 32, wherein the treatment 
candidate is one or more electrical impulses. 

36. The method according to claim 32, wherein said identifying 
identifies treatment candidates which achieve velocity restitution values consistent with 
histograms b and c of Figure 4 in "Dynamic Mechanism for Conduction Block in Heart 
Tissue" as therapeutic strategies for prevention of ventricular fibrillation. 

37. The method according to claim 32, wherein said identifying 
identifies treatment candidates which achieve potential duration restitution values 
consistent with histograms b and c of Figure 5 in "Dynamic Mechanism for Conduction 
Block in Heart Tissue" as therapeutic strategies for prevention of ventricular fibrillation. 
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38. The method according to claim 32, said identifying identifies 
treatment candidates which achieve cardiac memory values consistent with histograms b 
and c of Figure 6 in "Dynamic Mechanism for Conduction Block in Heart Tissue" as 
therapeutic strategies for prevention of ventricular fibrillation. 
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